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Abstract 

We study a heavy piston of mass M that separates finitely many ideal, 
unit mass gas particles moving in two or three dimensions. Neishtadt and 
Sinai previously determined a method for finding this system's averaged equa- 
tion and showed that its solutions oscillate periodically. Using averaging 
techniques, we prove that the actual motions of the piston converge in prob- 
ability to the predicted averaged behavior on the time scale M 1 / 2 when M 
tends to infinity while the total energy of the system is bounded and the 
number of gas particles is fixed. 
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1 Introduction 

Consider the following simple model of an adiabatic piston separating two gas 
containers: A massive piston of mass M 3> 1 divides a container in M. d into two 
halves. The piston has no internal degrees of freedom and can only move along 
one axis of the container. On either side of the piston there are a finite number of 
ideal, unit mass, point gas particles that interact with the walls of the container 
and with the piston via elastic collisions. When M = oo, the piston remains fixed 
in place, and each gas particle performs billiard motion at a constant energy in 
its sub- container. We make an ergodicity assumption on the behavior of the gas 
particles when the piston is fixed. Then we study the motions of the piston when 
the number of gas particles is fixed, the total energy of the system is bounded, but 
M is very large. 
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Heuristically, after some time, one expects the system to approach a steady 
state, where the energy of the system is equidistributed amongst the particles and 
the piston. However, even if we could show that the full system is ergodic, an 
abstract ergodic theorem says nothing about the time scale required to reach such 
a steady state. Because the piston will move much slower than a typical gas particle, 
it is natural to try to determine the intermediate behavior of the piston by averaging 
techniques. By averaging over the motion of the gas particles on a time scale chosen 
short enough that the piston is nearly fixed, but long enough that the ergodic 
behavior of individual gas particles is observable, we will show that the system 
does not approach the expected steady state on the time scale M 1//2 . Instead, the 
piston oscillates periodically, and there is no net energy transfer between the gas 
particles. 

This paper follows earlier work by Neishtadt and Sinai [NS04, Sin99j. They de- 
termined that for a wide variety of Hamiltonians for the gas particles, the averaged 
behavior of the piston is periodic oscillation, with the piston moving inside an ef- 
fective potential well whose shape depends on the initial position of the piston and 
the gas particles' Hamiltonians. They pointed out that an averaging theorem due 
to Anosov [Ano60[ ILM88j . proved for smooth systems, should extend to this case. 
This paper proves that Anosov's theorem extends to the particular gas particle 
Hamiltonian described above. Thus, if we examine the actual motions of the piston 
with respect to the slow time r = t/M 1 ^ 2 , then, as M — > oo, in probability (with 
respect to Liouville measure) most initial conditions give rise to orbits whose ac- 
tual motion is accurately described by the averaged behavior for < r < 1, i.e. for 
< t < M 1 / 2 . Gorelyshev and Neishtadt |GN06| and we |Wri06j have already 
proved that when d — 1, i.e. when the gas particles move on a line, the convergence 
of the actual motions to the averaged behavior is uniform over all initial conditions, 
with the size of the deviations being no larger than 0(M~ 1 ^ 2 ) on the time scale 
M-V2. 

The system under consideration in this paper is a simple model of an adiabatic 
piston. The general adiabatic piston problem [Cal63j, well-known from physics, 
consists of the following: An insulating piston separates two gas containers, and 
initially the piston is fixed in place, and the gas in each container is in a separate 
thermal equilibrium. At some time, the piston is no longer externally constrained 
and is free to move. One hopes to show that eventually the system will come to a 
full thermal equilibrium, where each gas has the same pressure and temperature. 
Whether the system will evolve to thermal equilibrium and the interim behavior 
of the piston are mechanical problems, not adequately described by thermodynam- 
ics |Gru99j . that have recently generated much interest within the physics and 
mathematics communities. One expects that the system will evolve in at least two 
stages. First, the system relaxes toward a mechanical equilibrium, where the pres- 
sures on either side of the piston are equal. In the second, much longer, stage, the 
piston drifts stochastically in the direction of the hotter gas, and the temperatures 
of the gases equilibrate. See for example |GPL03t IOL02] and the references therein. 
So far, rigorous results have been limited mainly to models where the effects of 
gas particles recolliding with the piston can be neglected, either by restricting to 
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extremely short time scales [LSC021 ICLS02] or to infinite gas containers [Che04j . 

A recent study involving some similar ideas by Chernov and Dolgopyat |CD0 6aj 
considered the motion inside a two-dimensional domain of a single heavy, large gas 
particle (a disk) of mass M ^> 1 and a single unit mass point particle. They as- 
sumed that for each fixed location of the heavy particle, the light particle moves 
inside a dispersing (Sinai) billiard domain. By averaging over the strongly hyper- 
bolic motions of the light particle, they showed that under an appropriate scaling 
of space and time the limiting process of the heavy particle's velocity is a (time- 
inhomogeneous) Brownian motion on a time scale 0(M 1 ^ 2 ). It is not clear whether 
a similar result holds for the piston problem, even for gas containers with good 
hyperbolic properties, such as the Bunimovich stadium. In such a container the 
motion of a gas particle when the piston is fixed is only nonuniformly hyperbolic 
because it can experience many collisions with the flat walls of the container im- 
mediately preceding and following a collision with the piston. 

The present work provides a weak law of large numbers, and it is an open 
problem to describe the sizes of the deviations for the piston problem [UD06bJ. 
Although our result does not yield concrete information on the sizes of the devia- 
tions, it is general in that it imposes very few conditions on the shape of the gas 
container. Most studies of billiard systems impose strict conditions on the shape 
of the boundary, generally involving the sign of the curvature and how the corners 
are put together. The proofs in this work require no such restrictions. In particu- 
lar, the gas container can have cusps as corners and need satisfy no hyperbolicity 
conditions. 

We begin in Section [2] by giving a physical description of our results. Precise 
assumptions and our main result, Theorem [T], are stated in Section [31 and a proof 
is presented in the following sections. 

2 Physical motivation for the results 

Before precisely stating our assumptions and results, we briefly review the physical 
motivations for our results and introduce some notation. 

Consider a massive, insulating piston of mass M that separates a gas container 
X> in R d , d = 2 or 3. See Fi gure [TJ Denote the location of the piston by Q, 
its velocity by dQ/dt = V, and its cross-sectional length (when d = 2, or area, 
when d = 3) by i. If Q is fixed, then the piston divides X> into two subdomains, 
T^i{Q) — T^i on the left and T> 2 (Q) = X> 2 on the right. By Ei we denote the total 
energy of the gas inside XV and by |X>j| we denote the area (when d = 2, or volume, 
when d = 3) of XV 

We are interested in the dynamics of the piston when the system's total energy 
is bounded and M — > oo. When M = oo, the piston remains fixed in place, and 
each energy Ei remains constant. When M is large but finite, MV 2 /2 is bounded, 
and so V = 0(M~ 1 ^ 2 ). It is natural to define 

e = M- 1 / 2 , W = -, 
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Figure 1: A gas container DcK 2 separated by a piston. 

so that W is of order 1 as e — > 0. This is equivalent to scaling time by e. 

If we let Pi denote the pressure of the gas inside T>i, then heuristically the 
dynamics of the piston should be governed by the following differential equation: 

§ -v, 

To find differential equations for the energies of the gases, note that in a short 
amount of time dt, the change in energy should come entirely from the work done 
on a gas, i.e. the force applied to the gas times the distance the piston has moved, 
because the piston is adiabatic. Thus, one expects that 

To obtain a closed system of differential equations, it is necessary to insert an 
expression for the pressures. Because the pressure of an ideal gas in d dimensions 
is proportional to the energy density, with the constant of proportionality 2/d, we 
choose to insert 

Pi 



dWA 



Later, we will make assumptions to justify this substitution. However, if we accept 
this definition of the pressure, and define the slow time 

t = at, 

we obtain the following ordinary differential equations for the four macroscopic 



4 



variables of the system: 
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(3) 



Neishtadt and Sinai |Sin99[ INS04j pointed out that the solutions of Equation (jSJ) 
have the piston moving according to an effective Hamiltonian. This can be seen as 
follows. Since 

d\V x {Q)\ . d\V 2 (Q)\ 



dQ 

dln{Ei)/dT = -(2/d)dln(\Vi(Q)\)/dr, and so 



dQ 
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Hence 



d 2 Q 2£ E^O) \V 1 {Q{0))\ 2/d 2£E 2 {0)\V 2 {Q{0))\ 2/d 



dr 2 d \V 1 (Q(r))\ 1+2/d d \V 2 (Q(r))\ 1+2/d ' 

and so (Q, W) behave as if they were the coordinates of a Hamiltonian system de- 
scribing a particle undergoing motion inside a potential well. The effective Hamil- 
tonian may be expressed as 



_ w2 , £i(0)|2MQ(0))| 



1 
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2/d 



+ 



E 2 (0)|X> 2 (Q(0))| 



i/d 



\v 2 {Q)\ 



2/d 



(4) 



The question is, do the solutions of Equation Q give an accurate description 
of the actual motions of the macroscopic variables when M tends to infinity? The 
main result of this paper is that, for an appropriately defined system, the answer 
to this question is affirmative for < t < M 1//2 , at least for most initial conditions 
of the microscopic variables. Observe that one should not expect the description 
to be accurate on time scales much longer than 0(M 1 ^ 2 ) = 0(e^ 1 ). The reason for 
this is that, presumably, there are corrections of size 0(e 2 ) in Equations (CQ) and 
(j2J) that we are neglecting. On the time scale e -1 , these errors roughly add up to 
no more than size 0(e^ 1 ■ e 2 = e), but on a longer time scale they should become 
significant. Such higher order corrections for the adiabatic piston were studied by 
Crosignani et al. [CDPS96J. 



3 Statement of the main result 
3.1 Description of the model 

We begin by describing the gas container. It is a compact, connected billiard 
domain Dcl d with a piecewise C 3 boundary, i.e. dT> consists of a finite number of 
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C 3 embedded hyper surf aces, possibly with boundary and a finite number of corner 
points. The container consists of a "tube," whose perpendicular cross-section V is 
the shape of the piston, connecting two disjoint regions. V C is a compact, 
connected domain whose boundary is piecewise C 3 . Then the "tube" is the region 
[0,l]xPcD swept out by the piston for < Q < 1, and [0, 1] x dP C dV. If 
d — 2, V is just a closed line segment, and the "tube" is a rectangle. If d — 3, V 
could be a circle, a square, a pentagon, etc. 
Our fundamental assumption is as follows: 

Main Assumption. For almost every Q e [0, 1] the billiard flow of a single particle 
on an energy surface in either of the two subdomains T>i(Q) is ergodic (with respect 
to the invariant Liouville measure). 

If d — 2, the domain could be the Bunimovich stadium [Bun79j . Another possible 
domain is indicated in Figure [U Polygonal domains satisfying our assumptions 
can also be constructed [Vor97j . Suitable domains in d = 3 dimensions can be 
constructed using a rectangular box with shallow spherical caps adjoined [BR98J. 
Note that we make no assumptions regarding the hyperbolicity of the billiard flow 
in the domain. 

The Hamiltonian system we consider consists of the massive piston of mass M 
located at position Q, as well as n\ + n 2 gas particles, n\ in T>i and n 2 in T> 2 . Here 
ni and n 2 are fixed positive integers. For convenience, the gas particles all have unit 
mass, though all that is important is that each gas particle has a fixed mass. We 
denote the positions of the gas particles in T>i by qij, 1 < j < rij. The gas particles 
are ideal point particles that interact with &D and the piston by hard core, elastic 
collisions. Although it has no effect on the dynamics we consider, for convenience 
we complete our description of the Hamiltonian dynamics by specifying that the 
piston makes elastic collisions with walls located at Q = 0, 1 that are only visible 
to the piston. We denote velocities by dQ/dt = V = eW and dqij/dt = v it j, and 
we set 

rii 

E id = vy2, /•:, ^ /•:,,, 

Our system has d(n± + n 2 ) + 1 degrees of freedom, and so its phase space is {2d{ni + 
n 2 ) + 2)-dimensional. 
We let 

h{z) = h = (Q, W, #1,1, E lt2 , • ■ ■ , E ltni ,E 2j i, E 2j2 , ■ • • , E 2t n 2 ), 

so that h is a function from our phase space to K n i+™2+2_ ^y e ft en abbreviate 
h = (Q, W, Eij, E 2 j), and we refer to h as consisting of the slow variables because 
these quantities are conserved when e = 0. We let h e (t, z) = h £ {t) denote the actual 
motions of these variables in time for a fixed value of e. Here z represents the initial 
condition in phase space, which we usually suppress in our notation. One should 
think of h £ (-) as being a random variable that takes initial conditions in phase space 
to paths (depending on the parameter t) in M. ni+Tl2+2 . 
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3.2 The averaged equation 

From the work of Neishtadt and Sinai [NS04], one can derive 

W 



d_ 



' Q ' 




w 


= H(h) := 







2E±l 2E 2 t 

d\Vi{Q)\ d\V 2 {Q)\ 
2WE 1 jt 



d\V 2 (Q)\ 



(5) 



as the averaged equation (with respect to the slow time r 
ables 



et) for the slow vari- 



Later, in Section 14.31 we will give another heuristic derivation of the aver- 
aged equation that is more suggestive of our proof. As in Section the solutions 
of Equation (jSJ) have (Q, W) behaving as if they were the coordinates of a Hamil- 
tonian system describing a particle undergoing motion inside a potential well. The 
effective Hamiltonian is given by Equation 
Let h(r, z) = h(r) be the solution of 



h(0) = h £ (0). 
Again, think of h(-) as being a random variable. 



3.3 The main result 

The solutions of the averaged equation approximate the motions of the slow vari- 
ables, h £ (t), on a time scale 0(1/ e) as e — > 0. Precisely, fix a compact set 
V C R n ^+ 2 such that h <E V => Q CC (0,1), W CC R, and E itj CC (0, oo) 
for each i and jo We will be mostly concerned with the dynamics when h E V. 
Define 

Qmin inf Qi Qmax SUp Q, 

hev heV 

I i 

E min = inf -W 2 + Ex + E 2 , E max = sup -W 2 + E 1 + E 2 . 
hev z /jgv * 

For a fixed value of e > 0, we only consider the dynamics on the invariant subset 
of phase space defined by 

M £ = {(Q, V, q itj , ViJ) e R2d(m+n 2) +2 : g e [o, i], % . e 

-Emm < ^"^ 2 + -^1 + -^2 < -Emax}- 

Let P £ denote the probability measure obtained by restricting the invariant Liouville 
measure to M. e . Define the stopping time 

T £ (z) =T £ = inf{r > : h(r) ^ V or h £ (r/e) £ V}. 



1 We have introduced this notation for convenience. For example, h € V Q CC (0, 1) means 
that there exists a compact set A C (0, 1) such that h E V Q G A, and similarly for the other 
variables. 
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Theorem 1. If V is a gas container in d = 2 or 3 dimensions satisfying the 
assumptions in Subsection \3.1\ above, then for each T > 0, 



sup \h £ (r/e) — h{r)\ —>■ in probability as e = M — > 0, 

0<t<TAT £ 



i. e. for each fixed 5 > 



Pe 



sup \h £ {r/e) — h{r)\ > 5 

,0<T<TAT e 




ase = M' 1/2 -> 0. 



It should be noted that the stopping time in the above result is not unduly 
restrictive. If the initial pressures of the two gasses are not too mismatched, then 
the solution to the averaged equation is a periodic orbit, with the effective potential 
well keeping the piston away from the walls. Thus, if the actual motions follow the 
averaged solution closely for < r < T A T £ , and the averaged solution stays in V, 
it follows that T £ > T. 

The techniques of this paper should immediately generalize to prove the ana- 
logue of Theorem [T] above in the nonphysical dimensions d > 3, although we do not 
pursue this here. 



4 Preparatory material concerning a two-dimensional 
gas container with only one gas particle on each 



Our results and techniques of proof are essentially independent of the dimension 
and the fixed number of gas particles on either side of the piston. Thus, we focus 
on the case when d = 2 and there is only one gas particle on either side. Later, 
in Section [6l we will indicate the simple modifications that generalize our proof to 
the general situation. For clarity, in this section and next, we denote q 1; i by qi, i>2,i 
by t>2, etc. We decompose the gas particle coordinates according to whether they 
are perpendicular to or parallel to the piston's face, for example q± = (gj 1 , gf ). See 
Figure El 

The Hamiltonian dynamics define a flow on our phase space. We denote this 
flow by z e (t,z) = z e (t), where z = z £ (0,z). One should think of z £ (-) as being a 
random variable that takes initial conditions in phase space to paths in phase space. 
Then h £ (t) = h(z £ (t)). By the change of coordinates W = V/e, we may identify all 
of the M. £ defined in Section [3] with the space 



and all of the P £ with the probability measure P on M. , which has the density 



side 



M = {(Q,W,q 1 ,v 1 ,q 2 ,v 2 ) GR 



Q e [0,1], q 1 eV 1 (Q), q 2 eV 2 (Q), 
E mm <-W* + E 1 + E 2 <E max }. 



dP = const dQdWdq 1 dq\dv l dv[dq 2 dq\dv 2 dv\. 
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Figure 2: A choice of coordinates on phase space. 

(Throughout this work we will use const to represent generic constants that are 
independent of e.) We will assume that these identifications have been made, so 
that we may consider z e (-) as a family of measure preserving flows on the same 
space that all preserve the same probability measure. We denote the components 
of z £ (t) by Q e (t), ql e (t), etc. 

The set {z G M. : q\ = Q = q 2 } has co-dimension two, and so \J t z e {t){qi = 
Q = q 2 } has co-dimension one, which shows that only a measure zero set of initial 
conditions will give rise to three particle collisions. We ignore this and other mea- 
sures zero events, such as gas particles hitting singularities of the billiard flow, in 
what follows. 

Now we present some background material, as well as some lemmas that will 
assist us in our proof of Theorem HJ We begin by studying the billiard flow of a gas 
particle when the piston is infinitely massive. Next we examine collisions between 
the gas particles and the piston when the piston has a large, but finite, mass. Then 
we present a heuristic derivation of the averaged equation that is suggestive of our 
proof. Finally we prove a lemma that allows us to disregard the possibility that 
a gas particle will move nearly parallel to the piston's face - a situation that is 
clearly bad for having the motions of the piston follow the solutions of the averaged 
equation. 

4.1 Billiard flows and maps in two dimensions 

In this section, we study the billiard flows of the gas particles when M = oo and 
the slow variables are held fixed at a specific value h G V. We will only study 
the motions of the left gas particle, as similar definitions and results hold for the 
motions of the right gas particle. Thus we wish to study the billiard flow of a point 
particle moving inside the domain T>i at a constant speed \f2E\. The results of 
this section that are stated without proof can be found in [CM06a] . 

Let TT>\ denote the tangent bundle to T>\. The billiard flow takes place in the 
three-dimensional space }A\ = M 1 = {(qi,Vi) G TT>\ : q\ G T>%, \vi\ = y/2E%}/ ~. 
Here the quotient means that when qi G <9T>i, we identify velocity vectors pointing 
outside of T>i with those pointing inside V>\ by reflecting through the tangent line 
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to &D% at qi, so that the angle of incidence with the unit normal vector to dVi 
equals the angle of reflection. Note that most of the quantities defined in this 
subsection depend on the fixed value of h. We will usually suppress this dependence, 
although, when necessary, we will indicate it by a subscript h. We denote the 
resulting flow by y(t,y) = y(t), where y(0,y) = y. As the billiard flow comes 
from a Hamiltonian system, it preserves Liouville measure restricted to the energy 
surface. We denote the resulting probability measure by \i. This measure has 
the density dfi = dqidvi/ {2ti^2Ei |£>i|). Here dq± represents area on M 2 , and dv\ 
represents length on = { u i e ^ 2 : l u il = \/2£a}- 

There is a standard cross-section to the billiard flow, the collision cross-section 



^ — {(<?i) v i) £ TT>i : qi G &Dx, \vi\ = \ / 2Ei}/ ~. It is customary to parameterize 
Q by {x = (r,(p) : r G &D%, <p G [— 7r/2, +7r/2]}, where r is arc length and <p 
represents the angle between the outgoing velocity vector and the inward pointing 
normal vector to &D 1 . It follows that Q may be realized as the disjoint union 
of a finite number of rectangles and cylinders. The cylinders correspond to fixed 
scatterers with smooth boundary placed inside the gas container. If F : Q O is the 
collision map, i.e. the return map to the collision cross-section, then F preserves the 
projected probability measure v, which has the density dv = cos ipdipdr/ (2 |<9X>i|). 
Here \&D\\ is the length of &V\. 

We suppose that the flow is ergodic, and so F is an invertible, ergodic mea- 
sure preserving transformation. Because &D± is piecewise C 3 , F is piecewise C 2 , 
although it does have discontinuities and unbounded derivatives near discontinu- 
ities corresponding to grazing collisions. Because of our assumptions on T>\, the 
free flight times and the curvature of &D\ are uniformly bounded. It follows that 
if x ^ dQ U then F is different iable at x, and 

11 ^ mi const 

DF(x) < — -, 6 

11 v ; " ~ cos^(Fx) v ; 

where <p(Fx) is the value of the p coordinate at the image of x. 

Following the ideas in Appendix |Aj we induce F on the subspace Cl of Q 
corresponding to collisions with the (immobile) piston. We denote the induced 
map by F and the induced measure by v. We parameterize Cl by {(r, ip) : < 
r < £, ip G [— 7r/2, +7r/2]}. As vVt = E/\dT>i\, it follows that v has the density 
dv = cos ipdipdr / (2£) . 

For x G define (x to be the free flight time, i.e. the time it takes the billiard 
particle traveling at speed \f2E[ to travel from x to Fx. If x ^ dVt U F~ 1 {dVl) 1 

Santalo's formula |San76[ IChe97j tells us that 
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If ( : Cl — > R is the free flight time between collisions with the piston, then it follows 
from Proposition [TU1 that 

The expected value of VjH when the left gas particle collides with the (immobile) 
piston is given by 



Eo \vH = E„\j2E~icos(p = / cos 2 <pd<p = v / 2£it- (10) 

2 J- n / 2 4 



We wish to compute lim^oot -1 f Q |2w4-(s)| 5 q ±^ = Qds, the time average of the 
change in momentum of the left gas particle when it collides with the piston. If this 
limit exists and is equal for almost every initial condition of the left gas particle, 
then it makes sense to define the pressure inside T>\ to be this quantity divided by 
I. Because the collisions are hard-core, we cannot directly apply Birkhoff 's Ergodic 
Theorem to compute this limit. However, we can compute this limit by using the 
map F. 

Lemma 2. // the billiard flow y(t) is ergodic, then for ^ — a.e. y G M 1 , 

1 I I / \l r , Exl 



1 f 

lim - / \v^(s)\6 q u s)=Q ds 

t— >oo I J 



2|Pi(Q)|' 



Proof. Because the billiard flow may be viewed as a suspension flow over the col- 
lision cross-section with ( as the height function, it suffices to show that the con- 
vergence takes place for v — a.e. x G Cl. For an initial condition x G Cl, define 
N t (x) = N t = # js G (0,t] : y(s,x) G Clj. By the Poincare Recurrence Theorem, 



N t — > oo as t — > oo, v — a.e. 
But 



i * i r* 

J2K\(F n x)< 7 / K(s)\5 qHs)=Q ds 

V < n=l 1 Jo 



N 1 Nt 

— sr^N.-l %/ x AT ^ I 1 I v ' 



and so the result follows from Birkhoff 's Ergodic Theorem and Equations (Q and 

(HDD. □ 

Corollary 3. // the billiard flow y(t) is ergodic, then for each 5 > 0, 



fi^yeM 1 : \vi(s)\8 q ± (s)=Q ds 



E x i 



> 5> — > as £ — > oo. 
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4.2 Analysis of collisions 

In this section, we return to studying our piston system when e > 0. We will 
examine what happens when a particle collides with the piston. For convenience, 
we will only examine in detail collisions between the piston and the left gas particle. 
Collisions with the right gas particle can be handled similarly. 

When the left gas particle collides with the piston, v ^ and V instantaneously 
change according to the laws of elastic collisions: 



1 — M 2M 
2 M - 1 





i 




l + M 





K1 




y- _ 



In our coordinates, this becomes 





1 


~e 2 -l 2e ' 




K1 


W+ 


l + e 2 


2e l-e\ 







Recalling that V\, W = 0(1), we find that to first order in e, 

vi+ = -v£- + 0(e), W + = W~ + 0(e) 



(12) 



Observe that a collision can only take place if v l 
-E\j2E max . Thus, either v\~ > or v\~ = 0(e). 
second order in e, it follows that 



> eW . In particular, v 1 > 
By expanding Equation ffTTl) to 



Et - E- x 
W + - W~ 



-2eW \vi\ +0(e 2 ) 
+2e\vH+0(e 2 ). 



(13) 



Note that it is immaterial whether we use the pre-collision or post-collision values 
of W and \v±\ on the right hand side of Equation ( fl3l) . because any ambiguity can 
be absorbed into the 0(e 2 ) term. 

It is convenient for us to define a "clean collision" between the piston and the 
left gas particle: 

Definition 1. The left gas particle experiences a clean collision with the piston if 
and only if > and v i + < —E\/2E max . 

In particular, after a clean collision, the left gas particle will escape from the pis- 
ton, i.e. the left gas particle will have to move into the region q± < before it can 
experience another collision with the piston. It follows that there exists a constant 
C\ > 0, which depends on the set V, such that for all e sufficiently small, so long 
as Q > Qmin and \v^\ > eC\ when q± 6 [Q m in,Q], then the left gas particle will 
experience only clean collisions with the piston, and the time between these colli- 
sions will be greater than 2Q min / (^2E max ). (Note that when we write expressions 
such as qi G [Qmin,Q], we implicitly mean that q± is positioned inside the "tube" 
discussed at the beginning of Section [3)) One can verify that C\ = 5\/2E max would 
work. 
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Similarly, we can define clean collisions between the right gas particle and the 
piston. We assume that C\ was chosen sufficiently large such that for all e suffi- 
ciently small, so long as Q < Q max and \v^\ > zC\ when e [Q, Q max \i then the 
right gas particle will experience only clean collisions with the piston. 

Now we define three more stopping times, which are functions of the initial 
conditions in phase space. 



T' £ =inf{r > : Q min < q^r/e) < Q e (r/e) < Q max and \v^ E (r/e)\ < C x e}, 
T e " =inf{r > : Q min < Q e {r/e) < g£ e (r/e) < Q max and \v^{t/s)\ < de}, 
f £ =T A T £ A T' £ A T' e ' 



Here we make use of Dirac delta functions. All integrals involving these delta 
functions may be replaced by sums. 

The following lemma is an immediate consequence of Equation (1131) and the 
above discussion: 

Lemma 4. // < t\ < t2 < T £ /e, the piston experiences 0((t2 — t\) V 1) colli- 
sions with gas particles in the time interval [tx,^], all of which are clean collisions. 
Furthermore, 



Here any ambiguities arising from collisions occurring at the limits of integration 
can be absorbed into the 0(e) term. 

4.3 Another heuristic derivation of the averaged equation 

The following heuristic derivation of Equation ([3]) when d = 2 was suggested 
in |Dol05j . Let At be a length of time long enough such that the piston expe- 
riences many collisions with the gas particles, but short enough such that the slow 
variables change very little, in this time interval. From each collision with the left 
gas particle, Equation ffl3|) states that W changes by an amount +2e \v±\ + 0(e 2 ), 
and from Equation ffTU]) the average change in W at these collisions should be ap- 
proximately ET\\j2E\j2 + 0(e 2 ). From Equation the frequency of these collisions 
is approximately \f2E\l/('K \T>\\). Arguing similarly for collisions with the other 
particle, we guess that 



Define H(z) by 



W 

+2\v^\S q ± =Q -2\v 2 L \S q ± =Q 
-2W vi 5 qi=Q 
+2W S q ^ =Q 




AW 



E 2 £ 



+ 0{e 2 ). 



At 



= e 



T>i(Q)\ 



— e 



v 2 (Q)\ 
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With r = et as the slow time, a reasonable guess for the averaged equation for W 
is 

dW E x l E 2 £ 

~dT~ \Vt{Q)\ ~ \v 2 (Q)\' 

Similar arguments for the other slow variables lead to the averaged equation (jSJ), 
and this explains why we used Pj = Eif \T>^\ for the pressure of a 2-dimensional gas 
in Section [3 

There is a similar heuristic derivation of the averaged equation in d > 2 dimen- 
sions. Compare the analogues of Equations © and (fTU]) in Subsection 16.21 



4.4 A priori estimate on the size of a set of bad initial 
conditions 

In this section, we give an a priori estimate on the size of a set of initial conditions 
that should not give rise to orbits for which sup 0<T<TATe \h e (r/e) — h(r)\ is small. 
In particular, when proving Theorem [Tj, it is convenient to focus on orbits that only 
contain clean collisions with the piston. Thus, we show that P{T £ < TAT £ } vanishes 
as e — ► 0. At first, this result may seem surprising, since P{T' e A T" = 0} = 0(e), 
and one would expect U t l z £ (— t){T' £ A T" = 0} to have a size of order 1. However, 
the rate at which orbits escape from {T' e A T" = 0} is very small, and so we can 
prove the following: 

Lemma 5. 

P{f £ <TAT £ } = 0(e). 

In some sense, this lemma states that the probability of having a gas particle 
move nearly parallel to the piston's face within the time interval [0, T/e], when 
one would expect the other gas particle to force the piston to move on a macro- 
scopic scale, vanishes as e — * 0. Thus, one can hope to control the occurrence of 
the "nondiffusive fluctuations" of the piston described in [CD06aJ on a time scale 

Proof. As the left and the right gas particles can be handled similarly, it suffices to 
show that P{T' e <T} = 0(e). Define 

23 e = {z G M : Qmin <qi<Q< Qmax and |^| < C x e}. 

Then {T' £ < T} C uJl £ z £ (-t)<B £ , and if 7 = Q min /V8E^, 

Me \ Me \ j Tie 

P |J Z e {-t)*8 e )=P\\J «e(*)»e = P » e U U((z e (*)8 e )\» e 
V t=0 / \ t=0 J \ t=0 



<P<B £ + Pl [j « e (fe7)[|j(%(*)»s)\» E 
V fe=0 t=o 

< P<B £ +(L + ij P ^U(^(t)» B )\» e j 
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Now PQ3 e = 0(e), so if we can show that P (ULofeW^A 9 ^) = °( £2 )' then [t 
will follow that P{T e ' < T} = 0(e). 

If z G UJLofeW^A^' ^ * s s t m true ^ na ^ l^i" = This is because [t^l 

changes by at most O(e) at the collisions, and if a collision forces |t>n > Cje, then 
the gas particle must escape to the region q± < before f j 1 can change again, and 
this will take time greater than 7. Furthermore, if z G UtLofeW^A^' then at 
least one of the following four possibilities must hold: 

• |?f - Qmin\ < 0(e), 

• IQ-Qminl < 0(e), 

• |Q - Qma^l < 0(e), 

• iQ-gf-l <0(e). 
It follows that P(UJ=o feW 



(9(e 2 ). For example, 



1 



{ ^ <0(e), 9j L 



const 



<o( £ )} 



dP 



{E mm <W 2 /2+vl/2+vl/2<E max } 
X / 1 



{Qe[o,iUie£>i,g2ez>2} 



{|i) ± |<0(£)}^^' r ^' y l ^1^2 ^2 



0{e 2 



□ 



5 Proof of the main result for two-dimensional 
gas containers with only one gas particle on 
each side 

As in Section HI we continue with the case when d = 2 and there is only one gas 
particle on either side of the piston. 

5.1 Main steps in the proof of convergence in probability 

By Lemma [5j it suffices to show that sup 0<T< j, £ \ h e (r/e) — h(r) \ — > in probability 
as e = M~ 1 / 2 — > 0. Several of the ideas in the steps below were inspired by 
a recent proof of Anosov's averaging theorem for smooth systems that is due to 
Dolgopyat |Dol05j . 
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Step 1: Reduction using Gronwall's Inequality. Observe that h{r) satisfies 
the integral equation 

h(r) - h(0) = [ H(h(a))da, 



while from Lemma HI 



h E (r/e) - h E (0) = 0(e) + e / H{z £ {s))ds 



/•t/s pt 

0{e)+e H{z e {s))- H{h E {s))ds+ H(h e {a/s))da 
Jo Jo 



for < r < T £ . Define 

e e (r) = e [ ' H(z £ (s))-H(h £ (s))ds. 
Jo 

It follows from Gronwall's Inequality that 

sup \h £ (r/e)-h(r)\< (o(e)+ sup \e £ (r)\) e Up ^) T . 

0<T<f E V 0<T<f e 



(14) 



Gronwall's Inequality is usually stated for continuous paths, but the standard proof 
(found in [SV85j ) still works for paths that are merely integrable, and \h e {r/e) — h{j) | 
is piecewise smooth. 

Step 2: Introduction of a time scale for ergodization. Let L(e) be a real 
valued function such that L{e) — > oo, but L{e) <^ loge^ 1 , as e — > 0. In Section l5\2l 
we will place precise restrictions on the growth rate of L(e). Think of L(e) as 
being a time scale that grows as e — * so that ergodization, i.e. the convergence 
along an orbit of a function's time average to a space average, can take place. 
However, L(e) doesn't grow too fast, so that on this time scale z £ (t) essentially 
stays on the submanifold {h = h £ (0)}, where we have our ergodicity assumption. 
Set tk, e = kL(e), so that 



eL(e) 



sup |e e (r)| < 0(eL(e)) +e ^ 



0<t<T £ 



k=0 



H{z £ {s)) - H{h £ {s))ds 



(15) 



Step 3: A splitting according to particles. Now H(z) — H(h(z)) divides into 
two pieces, each of which depends on only one gas particle when the piston is held 
fixed: 



H(z)-H(h(z)) 







2 K I ^qi=Q |Di 

-2W\vi\5 q , =Q + mm 




+ 



E 2 e 
\v 2 (Q)\ 



l-wM\+ 2W H\s q ^ Q ] 





- 2 \v. 
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We will only deal with the piece depending on the left gas particle, as the right 
particle can be handled similarly. Define 



G{z) = \vi\8 qi 



G{h) 



2|2MQ)| 



(16) 



Returning to Equation (115j) . we see that in order to prove Theorem (TJ it suffices to 
show that both 



E 

/c=0 



tfc+i, 



G{z £ {s)) - G{h £ {s))ds 



f -k. 



and 



eL(e) 



-1 



E 

/c=0 



tk+l, 



W £ (s)(G(z £ (s))-G(h £ (s)))ds 



tfe, 



converge to in probability as e — > 0. 



Step 4: A splitting for using the triangle inequality. Now we let 2& e (s) be 
the orbit of the e = Hamiltonian vector field satisfying Zk <£ (tk <e ) = 2 e (tfc e ). Set 
hk,e(t) — h(zk )S (t))- Observe that hk )£ (t) is independent of t. 

We emphasize that so long as < t < T £ /e, the times between collisions of a 
specific gas particle and piston are uniformly bounded greater than 0, as explained 
before Lemma HI It follows that, so long as tk+i, e < T £ /e, 

sup \h k , £ (t)-h £ (t)\ = 0(eL(e)). (17) 

*fe,E<*<tfe+l,e 

This is because the slow variables change by at most 0(e) at collisions, and dQ £ / dt = 
0(e). 
Also, 

f tk+ "w £ (s)(G(z £ (s)) -G(h e (s)))ds 

J *fe,e 

= 0(eL(e) 2 ) + W Ke {t Ke ) / G(z £ (s)) - G(h £ (s))ds, 

and so 




W E (s)(G(z E (s))-G(h £ (s)))ds 



< 0{eL(e)) +e const ^ 

k=0 



■k+l,e 



G(z £ (s))-G(h £ (s))ds 



tk,s 
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Thus, in order to prove Theorem [T] it suffices to show that 



eL(e) 

e 

k=0 



-1 



<£ ^ + \II k , £ \ + \III k>£ 



G(z £ (s))-G(h £ (s))ds^ 
converges to in probability as e — > 0, where 

h,s= [ k+l £ G(z £ (s))-G(z kte (s))ds, 
Ih,e = f k+1 ' e G(z kjE (s)) - G(h kj£ (s))ds, 

IIh,e = [ tk+le G(h k ,e(s)) ~ G(h e (s))d S . 

The term IIk, £ represents an "ergodicity term" that can be controlled by our 
assumptions on the ergodicity of the flow zo(t), while the terms I k)£ and IIIk,e 
represent "continuity terms" that can be controlled by controlling the drift of z £ (t) 
from z kt£ (t) for t k>£ < t < t k+1>£ . 

Step 5: Control of drift from the e = orbits. Now G is uniformly Lipschitz 
on the compact set V, and so it follows from Equation (fl7|) that III k , £ = 0(eL(e) 2 ). 



Thus, e J2 e k =o \IHk,e\ = 0{eL{e)) -> as e -> 0. 

Next, we show that for fixed S > 0, P I e J2k=o \^k,e\ — $ ) — * as £ ^ 0. 
For initial conditions z G M. and for integers k G [0, T/{eL[e)) — 1] define 



1 , 5 , . T, 



A k . £ = < z : — - |J fe ,| > — and < 



L(e) 1 2T - eL(e) 

A Zj£ = {k: z G • 

Think of these sets as describing "poor continuity" between solutions of the e = 
and the e > Hamiltonian vector fields. For example, roughly speaking, z G A kj£ 
if the orbit z £ {t) starting at z does not closely follow z k>£ (t) for t k>£ <t< t k+ i >£ . 

One can easily check that \I k ,e\ < ^i^{ £ )) f° r & ^ T £ /(eL(e)) — 1, and so it 
follows that 

Jk i 

£ \Ik,e\<- + 0(eL(e)#(A z , £ )). 

k=0 

Therefore it suffices to show that P (#(A Z , S ) > S (const eL^e))^ 1 ) — > as e — > 0. 
By Chebyshev's Inequality, we need only show that 

^ i 

E P (eL{e)#{A z J) = eL{e) ]T P(A, £ ) 

fc=0 



tends to with e. 

Observe that z £ (tk )£ )Ak, £ C Ao >E . In words, the initial conditions giving rise 
to orbits that are "bad" on the time interval [tk,£,tk+i, £ ], moved forward by time 
tk,e, are initial conditions giving rise to orbits which are "bad" on the time interval 
[to,e,ti,e\- Because the flow z e {-) preserves the measure, we find that 



T 

eL(s)- 



eL{e) P(A k>e ) < const P(Ao,e) 



k=0 

To estimate P(Ao t£ ), it is convenient to use a different probability measure, 
which is uniformly equivalent to P on the set {zeM: h(z) G V} D {T £ > eL(e)}. 
We denote this new probability measure by P* , where the / stands for "factor." If 
we choose coordinates on M. by using h and the billiard coordinates on the two gas 
particles, then P* is defined on Ai by dP* = dhd^d/j^, where dh represents the 
uniform measure on Vcl 4 , and the factor measure d[i\ represents the invariant 
billiard measure of the i th gas particle coordinates for a fixed value of the slow 
variables. One can verify that l{h{z)&v}dP < const dP$ , but that P^ is not invariant 
under the flow z £ (-) when e > 0. 

We abuse notation, and consider /a\ to be a measure on the left particle's initial 
billiard coordinates once h and the initial coordinates of the right gas particle are 
fixed. In this context, /x^ is simply the measure n from Subsection 14.11 Then 



P f (Ao, t 



< 



dh dn\ ■ n\ 



L(s) 



L{e) 



G{z £ {s)) - G(z Q (s))ds 



> ^ and eL{e) < f £ 



and we must show that the last term tends to with e. By the Bounded Conver- 
gence Theorem, it suffices to show that for almost every h G V and initial condition 
for the right gas particle, 



1 



L{e) 



G(z e {s)) - G(z (s))ds 



> — and eL{e) < f £ 



as e — > 0. 



Note that if G were a smooth function and z £ (-) were the flow of a smooth 
family of vector fields Z(z,e) that depended smoothly on e, then from Gronwall's 
Inequality, it would follow that sup < 4 < L(e) \z £ {t) - z (t)\ < 0(eL(e)e Lipiz)L(£) ). If 

this were the case, then L(s)~ l j Q L{£) G{z £ {s)) - G{z {s))ds = 0{eL{e)e Lip( - z ^), 

which would tend to with e. Thus, we need a Gronwall-type inequality for billiard 
flows. We obtain the appropriate estimates in Section 15.21 



Step 6: Use of ergodicity along fibers to control //, 



fc,e - 



be shown is that for fixed 5 > 0, P I e Ylk=o I-^mI — ^ 



-i 



All that remains to 
as e -> 0. 
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For initial conditions z G M and for integers k G [0,T/(eL(e)) — 1] define 



B. 



k,e 



B-. 




S T 
Ilk e I > and k < - e . . 
,e| 2T ~ eL(e) 



- 1 



Think of these sets as describing "bad ergodization." For example, roughly speak- 
ing, z G Bk, £ if the orbit z £ (t) starting at z spends the time between tk, e and tk+i, £ 
in a region of phase space where the function G(-) is "poorly ergodized" on the 
time scale L(e) by the flow zo(t) (as measured by the parameter 5/2T). Note that 
G(z) — \Vi\ 5 q ± = Q is not really a function, but that we may still speak of the con- 
vergence of t _1 J** G(zo(s))ds as t — > oo. As we showed in Lemma [2], the limit is 
G{Hq) for almost every initial condition. 

Proceeding as in Step 5 above, we find that it suffices to show that for almost 
every h G V, 



G(z (s))ds - G(h (0)) 



> — )> — > as t — > oo. 
~ 2T 



But this is simply a question of examining billiard flows, and it follows immediately 
from Corollary [3] and our Main Assumption. 



5.2 A Gronwall-type inequality for billiards 

We begin by presenting a general version of Gronwall's Inequality for billiard maps. 
Then we will show how these results imply the convergence required in Equa- 
tion (1181). 



5.2.1 Some inequalities for the collision map 

In this section, we consider the value of the slow variables to be fixed at ho G V. 
We will use the notation and results presented in Section H~T| but because the value 
of the slow variables is fixed, we will omit it in our notation. 

Let p, 7, and A satisfy 0<p^C7^Cl^cA<oo. Eventually, these quantities 
will be chosen to depend explicitly on e, but for now they are fixed. 

Recall that the phase space Q for the collision map F is a finite union of disjoint 
rectangles and cylinders. Let d be the Euclidean metric on connected components 
of Q. If x and x' belong to different components, then we set d(x,x') = oo. The 
invariant measure v satisfies v < const ■ (Lebesgue measure). For A C Q and a > 0, 
let Af a (A) = {x G Q : d(x, A) < a} be the a-neighborhood of A. 

For x G f2 let Xk(x) = x^ = F k x, k > 0, be its forward orbit. Suppose x ^ C 7i a, 
where 

c 7 ,a = (uL F- fc A/;(«9fi)) |J(uto^W 7 (F-W 7 (fln))). 

Thus for < k < A, is well defined, and from Equation (jH]) it satisfies 

d<y,£fc)<7 =^ d(Fx',x k+1 ) <^^-d(x',x k ). (19) 
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Next, we consider any p-pseudo-orbit x' k obtained from x by adding on an error 
of size < p at each application of the map, i.e. d(x' ,xo) < p, and for k > 1, 
d(x' k , Fx' k _]) < p. Provided d(xj, x'A < 7 for each j < k, it follows that 

. . ^-^ /const V /const \ k 

d(x k ,x k ) < p 2_l I < const pi 1 . (20) 

j=o V 7 / V 7 / 

In particular, if p, 7, and A were chosen such that 

/ const \ A ,„„ , 

const pi J < 7, (21) 

then Equation (I2"0"j) will hold for each k < X. We assume that Equation (|2~T|) is true. 
Then we can also control the differences in elapsed flight times using Equation ([7]): 

k 



_ , , const p / const \ . , 

\(x k -(x' k \< 'I . (22) 

7 V 7 / 



It remains to estimate the size z/C 7; a of the set of x for which the above estimates 
do not hold. Using Lemma M below, 

vC^ x < {\ + l)(vM J (d£l) + vtf J (F- 1 M J (da))) < 0(A( 7 + 7 1/3 )) = C(A 7 1/3 ). (23) 

Lemma 6. As 7 — > 0, 

i/JV 7 (F-W 7 (0fi)) = C(7 1/3 ). 

This estimate is not necessarily the best possible. For example, for dispersing 
billiard tables, where the curvature of the boundary is positive, one can show that 
vN^F -1 N^dVl)) = 0(7). However, the estimate in Lemma El is general and 
sufficient for our needs. 

Proof. First, we note that it is equivalent to estimate vJ\fy(FAfy(d£l)) , as F has the 
measure-preserving involution X(r, <p) = (r, —if), i.e. F~ x = IofoI |CM 06bj . 

Fix a G (0, 1/2), and cover J\fy(d£l) with C(7 _1 ) starlike sets, each of diameter 
no greater than 0(7). For example, these sets could be squares of side length 7. 
Enumerate the sets as {AA-. Set G = {i : FAi n J\tya(dfl) = 0}. 

li i E G, F\a { is a diffeomorphism satisfying HDF^JI < 0(j~ a ). See Equa- 
tion ©. Thus diameter (FA) < 0(7 1_a ), and so diameter (A/^FA;)) < C(7 1_a ). 
Hence vJ\f 7 (FAi) < O^ 1 " )), and z/A/^U^gFA) < C( 7 1 - 2a ). 

If z ^ £/, Ai (7 F~ 1 (A/" 7 «(9fi)) 7^ 0. Thus Ai might be cut into many pieces 
by F _1 (<9f2), but each of these pieces must be mapped near dVL. In fact, FAi C 
Afo^dQ). This is because outside F- l (Mf>(dSl)), \\DF\\ < C(7~ Q ), and so 
points in FAi are no more than a distance 0(7/7°) away from A/" 7 a(<9f2), and 
7 < 7 1 " < 7 a . It follows that Af^FAA C 7V 0(7Q) (9fi), and uU {^){d9) = 0(^ a ). 

Thus vM 1 {F~ l M 1 {dVt)) = (9(7 1_2q: + 7"), and we obtain the lemma by taking 
a = 1/3. 

□ 
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5.2.2 Application to a perturbed billiard flow 

Returning to the end of Step 5 in Section 15.11 let the initial conditions of the slow 
variables be fixed at ho = (Qo, Wo, i?i,o> E 2i o) G V throughout the remainder of this 
section. We can assume that the billiard dynamics of the left gas particle in T>i(Q ) 
are ergodic. Also, fix a particular value of the initial conditions for the right gas 
particle for the remainder of this section. Then z E (t) and T £ may be thought of 
as random variables depending on the left gas particle's initial conditions y G Ai l . 
Now if h e {t) = (Q e (t), W e (t), Ei tS (t), E2 >e {t)) denotes the actual motions of the slow 
variables when e > 0, it follows from Equation (fl7|) that, provided eL{e) < T e , 

sup \ho-h £ {t)\=0{eL{e)). (24) 

0<t<L(e) 



Furthermore, we only need to show that 
piyeM 1 : 



1 /L '~ G(z £ (s)) - G(z (s))ds 



m Jo 



> — and eL{s) < T £ } - 



(25) 

as e — > 0, where G is defined in Equation (fl6l) . 

For definiteness, we take the following quantities from Subsection 15.2.11 to de- 
pend on e as follows: 

L(e) =L = log log-, 

70) = 7 = e~ L , 

AW = A = ^, (26) 
E„(, 

, . eL 
p{E) = p = const — . 

7 

The constant in the choice of p and p's dependence on e will be explained in the 
proof of Lemma [U which is at the end of this subsection. The other choices may be 
explained as follows. We wish to use continuity estimates for the billiard map to 
produce continuity estimates for the flow on the time scale L. As the divergence of 
orbits should be exponentially fast, we choose L to grow sub logarithmically in e^ 1 . 
Since from Equation (jSJ) the expected flight time between collisions with &Di(Qo) 
when e = is E„( = tt \T> 1 (Q )\ /(y/2Ei$ \dV 1 (Q )\), we expect to see roughly 
A/2 collisions on this time scale. Considering A collisions gives us some margin for 
error. Furthermore, we will want orbits to keep a certain distance, 7, away from 
the billiard discontinuities. 7 — > as e — ► 0, but 7 is very large compared to the 
possible drift 0(eL) of the slow variables on the time scale L. In fact, for each 
C,m,n > 0, 

- = C(£e consti2 ) ^0 as^^ 0. (27) 

T \l J 

Let X : M. 1 — > Q be the map taking y £ Ai 1 to x = X(y) G Q, the location 
of the billiard orbit of y in the collision cross-section that corresponds to the most 
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recent time in the past that the orbit was in the collision cross-section. We consider 
the set of initial conditions 

s e = x-^v^a) H x_1 { x G n : ^ pkx ) > L 

Now from Equations (I2"3"j) and (1251) . vC^^ — > as e — > 0. Furthermore, by the ergod- 
icity ofF,v{xeQ: Eto C(^) < = z/ {x G : A" 1 £to C(^) < «/2} 
as e — > 0. But because the free flight time is bounded above, /iX™ 1 < const -u, and 
so \xE e — > 1 as e — > 0. Hence, the convergence in Equation (1251) and the conclusion 
of the proof in Section 15.11 follow from the lemma below and Equation (J2j 

Lemma 7 (Analysis of deviations along good orbits). As e —>■ 0, 

-L 



sup 

y&£ s n{eL<f e } 



A 

0(L- 1 ) -> 0. 



Proof. Fix a particular value of y E £ £ fl jeL < T e |. For convenience, suppose that 

y = X(y) = x G Q. Let yo{t) denote the time evolution of the billiard coordinates 
for the left gas particle when e = 0. Then there is some N < A such that the 
orbit Xk = F k x = (r^, (p^) for < k < N corresponds to all of the instances 
(in order) when yo(t) enters the collision cross-section Q = Qh corresponding to 
collisions with &Dx(Qq) for < t < L. We write Qh to emphasize that in this 
subsection we are only considering the collision cross-section corresponding to the 
billiard dynamics in the domain T>i(Qq) at the energy level -Ei,o- In particular, F 
will always refer to the return map on Qh - 

Also, define an increasing sequence of times corresponding to the actual times 
yo(t) enters the collision cross-section, i.e. 

* = 0, 

tk = tfe-i + C^Jfc-i for k > 0. 
Then Xk = Vo(tk)- Furthermore, define inductively 

Ni = inf {k > : tk corresponds to a collision with the piston} , 
Nj = inf {k > Nj~i : tk corresponds to a collision with the piston} . 

Next, let y e {t) denote the time evolution of the billiard coordinates for the left 
gas particle when e > 0. We will construct a pseudo-orbit x' ke = (r' k ^,(p k ^ £ ) of 
points in Qh that essentially track the collisions (in order) of the left gas particle 
with the boundary under the dynamics of y e {t) for < t < L. 

First, define an increasing sequence of times t' ke corresponding to the actual 
times y e (t) experiences a collision with the boundary of the gas container or the 
moving piston. Define 

N' e = sup {k > : t' k £ < L) , 

N[ e = inf > : t' k e corresponds to a collision with the piston} , 

N'j E = inf {A; > iV'-_ ls : t' ke corresponds to a collision with the piston} . 



23 



X 



k.e 



Because L < T £ (y)/e, we know that as long as N' j+1>e < N' £ , then iVj +1)£ - N' j>£ > 2. 
See the discussion in Subsection 14.21 Then we define x' k £ G f^ by 

'y £ (tU if^KJ, 
F-'x' k+h£ if ke{N>J. 

Lemma 8. Provided e is sufficiently small, the following hold for each k G [0, N A 
N' e ). Furthermore, the requisite smallness of e and the sizes of the constants in these 

estimates may be chosen independent of the initial condition y G £ £ H jeL < T e 

and of k: 



(a) x't 



is we 



11 defined. In particular, if k {A^ e } ; y £ (t' k£ ) corresponds to 



collision point on dT>i(Q ), and not to a collision point on a piece of dT) to 
the right of Qq. 

(b) Ifk>0andk(£ {N' j£ }, then x' k£ = Fx' k _ lE . 

(c) If k > and k G {iV' e | ; then d(x' k e , Fx' k _ l £ ) < p and the ip coordinate of 
y £ {t' ke ) satisfies <p(y £ (t kjE )) = <f' k)E + 0(e). 

(d) d(xk,x ke ) < const p(const/^) k . 

(e) k = iVj j£ if and only if k = Nj . 



(f) lfk>0,t>-t> k 



k-l,e 



t k - t k -! + 0(p (const/j) 



We defer the proof of Lemma [8] until the end of this subsection. Assuming that 
e is sufficiently small for the conclusions of Lemma [8] to be valid, we continue with 
the proof of Lemma [7J 

Set M = N A N £ - 1. Note that M < A ~ L. From (f) in Lemma E and 
Equations fl2T)|) and f[2"Tj) . we see that 



as e — > 0. 



/ I I ' ' I /n ( cons ^ A \ 

t M - t Mj£ \ < \H,e ~ h-l,e ~ {tk - tk-l)\ = O I p j — 1 

fe=i V 7 / 



Because the flight times t' k£ — t' k _ l£ and t k — t k -\ are uniformly bounded above, 
it follows from the definitions of iV and N' £ that tM, t' M £ > L — const. But from 
Subsection I4.2[ the time between the collisions of the left gas particle with the 
piston are uniformly bounded away from zero. Using (c) and Equation (|24|) . it 
follows that 



G{z £ {s)) - G{z {s))ds 
OiL- 1 ) + 

keiNf.NjKM} 
ke{Ny.Nj<M} 

OiL' 1 ) + 0(eL 2 ) + ^2E^ 

ke{Ny.N 3 <M} 



2E lfl cosy^fc - J2E h£ (t' k£ ) cosfaL + O(e)) 



2E lfi costp k - y/2E 1<0 costp' k£ + 0{eL) 



COSy? fc - cos^ fc 
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But using (d), 

M 

|cos<£>fc — cosy4 je | < (9(p(const/7) fc ) = (9(p(const/7) A ). 

kE{Nj-.Nj<M} k=0 

Since eL 2 = (9(p(const/7) A ), this finishes the proof of LemmaEl 

□ 

Proof of Lemma The proof is by induction. We take e to be so small that Equa- 
tion (|2ip is satisfied. This is possible by Equation (1271) . 

It is trivial to verify (a)-(f) for k = 0. So let < i < N A A^, and suppose that 
(a)-(f) have been verified for all k < I. We have three cases to consider: 

Case 1: / - 1 and / £ {iVj J: 

In this case, verifying (a)-(f) for k = I is a relatively straightforward application 
of the machinery developed in Subsection I5.2.1[ because for t' l _ l£ <t<t' le , y E (t) 
traces out the billiard orbit between x' l _ l£ and x' le corresponding to free flight in 
the domain T> 1 (Q ). We make only two remarks. 

First, as long as e is sufficiently small, it really is true that x' le = y £ (t' le ) corre- 
sponds to a true collision point on dT>i(Q ). Indeed, if this were not the case, then 
it must be that Q £ (t' le ) > Qo, and y £ {t' l£ ) would have to correspond to a collision 
with the side of the "tube" to the right of Qo- But then x" £ = Fx' l _ l£ G O^ 
would correspond to a collision with an immobile piston at Q and would satisfy 
d(xk,x'l, e ) < const p( const /7) fc < const p(const/7) A = 0(7), using Equations fl2"Ul) 
and (12"T|) . But x^ ^ A/" 7 (<9i\ ), and so it follows that when the trajectory of y £ (t) 
crosses the plane {Q = Qo}, it is at least a distance ~ 7 away from the boundary of 
the face of the piston, and its velocity vector is pointed no closer than ~ 7 to being 
parallel to the piston's face. As Q £ {t' le ) — Q = (D(sL) = 0(7), and it is geometri- 
cally impossible (for small e) to construct a right triangle whose sides si, S2 satisfy 
>~ 7> I s 2 1 < 0(eL), with the measure of the acute angle adjacent to Si being 
greater than ~ 7, we have a contradiction. After crossing the plane {Q = Qo}, 
y £ (t) must experience its next collision with the face of the piston, which violates 
the fact that / £ {N' jj£ }. 

Second, t' le — t' l _ l£ = C x 'i-i. £ + 0(eL), because t>i i£ = t>i i0 + 0(eL). See 
Equation ([24}. From Equation [22l |C^/-i - C^-i, e | < 0({pll) ( const /t)' -1 )- As 
ti — = (xi-i and sL = 0((p/-y) (const/7)' -1 ), we obtain (f). 

Case 2: There exists i such that / = N' i£ : 

For definiteness, we suppose that Q e (t; £ ) > Qo, so that the left gas particle collides 
with the piston to the right of Qo- The case when Q £ (t' le ) < Q can be handled 
similarly. 
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We know that Xi, x\+\ $L Afj(d{l ho ) U Af y (F W 7 (9fi/ l0 )). Using the induc- 
tive hypothesis and Equation (I2(jp . we can define 

r " - TP r i r " - p2 r / 

x l,e ~ r x l~l,e' x l+l,e ~ r x Z-l,e> 

and d(xi,Xi e ) < const p(const/7)', d(xi + \, x" +1 £ ) < const p(const/7)' +1 . In particu- 
lar, x'( £ and x'{ +1 E are both a distance ~ 7 away from dQh - Furthermore, when the 
left gas particle collides with the moving piston, it follows from Equation (j!2p that 
the difference between its angle of incidence and its angle of reflection is 0(e). Re- 
ferring to Figure [31 this means that ip\ £ = Lp'l £ + 0{e). Geometric arguments similar 
to the one given in Case 1 above show that the ^-trajectory of the left gas particle 
has precisely one collision with the piston and no other collisions with the sides of 
the gas container when the gas particle traverses the region Q < Q < Q £ (t' l£ ). Note 
that x\ e was defined to be the point in the collision cross-section Qh corresponding 
to the return of the y e -trajectory into the region Q < Q . See Figure [3j From this 
figure, it is also evident that d{r[ £ , r" e ) < 0(sL/j). Thus d(x' l ' e , x\ £ ) = 0(eL/ r y), 
and this explains the choice of p(e) in Equation (T2"6"|) . 

From the above discussion and the machinery of Subsection I5.2.1[ (a)-(e) now 
follow readily for both k = I and k — I + 1. Furthermore, property (f) follows in 
much the same manner as it did in Case 1 above. However, one should note that 
t[ £ -t[_^ = Cx'^ + O^ + O^eLh) andt' l+h£ -t[ £ = C^ + 0(eL)+0(eL/ 7 ), 
because of the extra distance 0{eLH) that the gas particle travels to the right of 
Qo- But eL/j = (9((p/7) (const/7) ~ )> an d so property (f) follows. 

Case 3: There exists i such that I — 1 — N' ie : 

As mentioned above, the inductive step in this case follows immediately from our 
analysis in Case 2. 

□ 

6 Generalization to a full proof of Theorem d 

It remains to generalize the proof in Sections H] and [5] to the cases when n\,n2>l 
and d = 3. 

6.1 Multiple gas particles on each side of the piston 

When d = 2, but rii,n 2 > 1, only minor modifications are necessary to generalize 
the proof above. As in Subsection 14.21 one defines a stopping time T £ satisfying 

P |t £ < T A T £ | = 0{e) such that for < t < T e /e, gas particles will only experi- 
ence clean collisions with the piston. 
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z>i(Qo) 



r— coordinate 




7 /2 



7 /2 



0(eL/ 7 ) 



0{eL) 



Figure 3: An analysis of the divergences of orbits when e > and the left gas particle 
collides with the moving piston to the right of Qq. Note that the dimensions are 
distorted for visual clarity, but that eL and eL/^y are both 0(7) as e — > 0. 
Furthermore, G (-tt/2+ 7 /2, tt/2- 7 /2) and = ^ e +0(e), and so r[ >£ = 
G{eL/ r y). In particular, the y e -trajectory of the left gas particle has precisely one 
collision with the piston and no other collisions with the sides of the gas container 
when the gas particle traverses the region Q < Q < Q £ {t' l £ ) 
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Next, define H(z) by 



H{z) = 







W 








=Q ~ 






-2W 


K 






+2W 







It follows that for < t < f e /e, h e (t) - h £ {0) = 0(e)+e jj H(z £ (s))ds. From here, 
the rest of the proof follows the same steps made in Subsection 15. 11 We note that 
at Step 3, we find that H(z) — H(h(z)) divides into n\ + n 2 pieces, each of which 
depends on only one gas particle when the piston is held fixed. 



6.2 Three dimensions 

The proof of Theorem [I] in d = 3 dimensions is essentially the same as the proof in 
two dimensions given above. The principal differences are due to differences in the 
geometry of billiards. We indicate the necessary modifications. 

In analogy with Section I4.1[ we briefly summarize the necessary facts for the 
billiard flows of the gas particles when M = oo and the slow variables are held 
fixed at a specific value h G V. As before, we will only consider the motions of one 
gas particle moving in T>\. Thus we consider the billiard flow of a point particle 
moving inside the domain T>\ at a constant speed \J2E\. Unless otherwise noted, 
we use the notation from Section 14.11 

The billiard flow takes place in the five- dimensional space M 1 = G 
TV\ : qi G T>i, \v\\ = \f2E\}/ ~. Here the quotient means that when qi G dVi, 
we identify velocity vectors pointing outside of T>\ with those pointing inside D\ 
by reflecting orthogonally through the tangent plane to &Di at q\. The billiard 
flow preserves Liouville measure restricted to the energy surface. This measure has 
the density d\i = dq\dv\/ (RiiEi \T>i\). Here dqi represents volume on M 3 , and dv\ 
represents area on S^^- = |t>i G IR 3 : \v\\ = ^2Ei}. 

The collision cross-section Q = {(q\,vi) G TT>\ : q\ G dV 1 , \v\\ = \f2E{}/ ~ is 
properly thought of as a fiber bundle, whose base consists of the smooth pieces of 
dT>i and whose fibers are the set of outgoing velocity vectors at q% G &D\. This and 
other facts about higher-dimensional billiards, with emphasis on the dispersing case, 
can be found in [BCST03"] . For our purposes, Q can be parameterized as follows. We 
decompose dT>i into a finite union U^Tj of pieces, each of which is diffeomorphic via 
coordinates r to a compact, connected subset of M 2 with a piecewise C 3 boundary. 
The Tj are nonoverlapping, except possibly on their boundaries. Next, if (qi,vi) G 
Q and v\ is the outward going velocity vector, let v = v\j \v\\. Then Q can be 
parameterized by {x = (r, £>)}. It follows that Q it is diffeomorphic to UjTj x 
S 2+ , where S 2+ is the upper unit hemisphere, and by dQ we mean the subset 
diffeomorphic to (U^ x S 2+ ) [JiUjTj x dS 2+ ). If x G fi, we let <p G [0,vr/2] 
represent the angle between the outgoing velocity vector and the inward pointing 
normal vector n to &D\, i.e. costp = (v, n). Note that we no longer allow (p to take 
on negative values. The return map F : Q O preserves the projected probability 
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measure v, which has the density dv = cos (p dv dr/(7r \ dT>i\). Here \&D\ | is the area 

F is an invertible, measure preserving transformation that is piecewise C 2 . Be- 
cause of our assumptions on T>\, the free flight times and the curvature of &D\ are 
uniformly bounded. The bound on given in Equation ([6]) is still true. A 

proof of this fact for general three-dimensional billiard tables with finite horizon 
does not seem to have made it into the literature, although see [BCST03] for the 
case of dispersing billiards. For completeness, we provide a sketch of a proof for 
general billiard tables in Appendix IBl 

We suppose that the billiard flow is ergodic, so that F is ergodic. Again, we 
induce F on the subspace Cl of fl corresponding to collisions with the (immobile) 
piston to obtain the induced map F : Cl O that preserves the induced measure v. 

The free flight time ( : Q — > K again satisfies the derivative bound given in 
Equation ([7j). The generalized Santalo's formula [Che97| yields 

EuC 



4|£>! 




Kl 





If C : Cl — > R is the free flight time between collisions with the piston, then it follows 
from Proposition [TU1 that 

\vi \ £ 

The expected value of mn when the left gas particle collides with the (immobile) 
piston is given by 

\vH = F^^plEx cosy? = ff cos 2 ip dv x = ^/2Fh-. 

TT JJ S 2+ 3 

As a consequence, we obtain 
Lemma 9. For fi — a.e. y £ A4 l , 



Compare the proof of Lemma [21 

With these differences in mind, the rest of the proof of Theorem [1] when d = 3 
proceeds in the same manner as indicated in Sections l4f l5l and 16. II above. The only 
notable difference occurs in the proof of the Gronwall-type inequality for billiards. 
Due to dimensional considerations, if one follows the proof of Lemma [6] for a three- 
dimensional billiard table, one finds that vM^{F~ l M^{dVl)) = 0(^ l ~ ia + 7°). The 
optimal value of a is 1/5, and so vN' 1 (F~ l N' 1 (dVl)) = O^ 1 ^) as 7 — * 0. Hence 
vC^^x = (9(A7 1 / 5 ), which is a slightly worse estimate than the one in Equation (123]) . 
However, it is still sufficient for all of the arguments in Section 15.2.2] and this 
finishes the proof. 
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A Inducing maps on subspaces 



Here we present some well-known facts on inducing measure preserving transforma- 
tions on subspaces. Let F : (fi, 23, v) O be an invertible, ergodic, measure preserv- 
ing transformation of the probability space Q endowed with the cr-algebra 23 and the 
probability measure v. Let Q G 23 satisfy < uCl < 1. Define R : Cl —>■ N to be the 
first return time to Cl, i.e. Ru = inf{n G N : F n u G Cl}. Then if v := v{- D f2)/i/£] 
and 23 := {F n : F G 23}, F : (fi,»,i>) O defined by Feu = F^u; is also 
an invertible, ergodic, measure preserving transformation |Pet83j . Furthermore 
EpR = J^Rdu = (z/fi)- 1 . 

This last fact is a consequence of the following proposition: 



Proposition 10. If ( : Q — > 1R> is in L 1 {v), then ( = J2 n =o C ° E n is in L^iy), 
and 

v\l 

Proof. 

» R— 1 » K—l oo „ k—l 

uQ / V C o F n dv = / V C o F n dv = V / V C o F n dv 

J n n=o J & n=0 k=1 Jnn{R=k} n=0 

oo k—l „ „ 

= )]V / Cdu= (dv, 

k=l n=0 JF"(Qn{R=k}) Jn 



because {F n (Cl fl {R = k}) : < n < k < oo} is a partition of fl 



□ 



B Derivative bounds for the billiard map in three 
dimensions 

Returning to Section 16. 2\ we need to show that for a billiard table T>\ C M 3 with 
a piecewise C 3 boundary and the free flight time uniformly bounded above, the 
billiard map F satisfies the following: If xq ^ dQ U F~ 1 (dfl), then 

ii ^ mi const 
\\DF(x )\\ < 



cosip(Fxo) 



Fix xo = (r ,vo) G f2, and let x\ = (r\,V\) = Fx®. Let S be the plane that 
perpendicularly bisects the straight line between r and r 1; and let r\j% denote 
the point of intersection. We consider £ as a "transparent" wall, so that in a 
neighborhood of xq, we can write F = F2 o F\. Here, F\ is like a billiard map in 
that it takes points (i.e. directed velocity vectors with a base) near xq to points with 
a base on £ and a direction pointing near r\. (Ft would be a billiard map if we 
reflected the image velocity vectors orthogonally through E.) F 2 is a billiard map 
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that takes points in the image of F\ and maps them near x\. Let x\ji = FiXq = 
F^W Then \\DF(x )\\ < ll^i(^o)|| ||DF 2 (z 1/2 )||. 

It is easy to verify that ||DFi(xo)|| < const, with the constant depending only on 
the curvature of &D\ at tq. In other words, the constant may be chosen independent 
oixo- Similarly, ||DF 2 _1 (a;i) || < const. Because billiard maps preserve a probability 
measure with a density proportional to cosy?, detZXF 2 ~ (xi) = cos^i/ cos(/?i/ 2 = 
cos^i. As Q is 4-dimensional, it follows from Cramer's Rule for the inversion of 
linear transformations that 

H np/ vii . const ||DF 2 -1 (xi) II 3 const 

\\DF-2 Xi/2 < u r- — — — < 

II V^ll — j„j-np-i/™ "\ — 



detDF 2 L (xi) coscpi' 



and we are done. 
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